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ABSTRACT 

The  OPBL  (Oceanic  Planetary  Boundary  Layer)  model  is  used  to  determine  the  rel- 
ative effects  of  the  seasonally  varying  surface  fluxes  of  heat,  salt  and  momentum  on  the 
seasonal  thermocline  and  mixed  layer  in  the  Mediterranean  Sea.  Three  combinations 
of  surface  forcing  are  tested  to  find  their  relative  contributions  to  mixing.  In  the  first 
case  free  convection  forcing  and,  in  the  second  case  wind  forced  convection  and 
buoyancy  flux  due  only  to  heat  are  considered  to  determine  their  individual  effects  on 
the  mixed  layer  deepening.  Then  in  the  third  case  both  free  and  wind  forced  convections 
are  applied  to  the  model  in  order  to  determine  their  combined  effect  on  the  mixed  layer 
depth. 

As  a  result,  evaporation  is  shown  to  have  a  significant  impact  on  the  mixed  layer 
deepening  and  on  the  temperature  and  salinity  structure  of  the  seasonal  pycnocline.  It 
has  been  shown  that  deepening  due  to  evaporation  exceeds  the  deepening  due  to  wind 
stress  between  October  and  January.  In  the  late  winter  the  strong,  cold,  dry  winds  from 
North  to  Northwest  (Mistral,  Tramontane)  causes  forced  convection  which  exceeds  the 
deepening  due  to  evaporation  at  this  time  of  the  year. 
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I.      INTRODUCTION 

A.  PURPOSE  OF  THE  STUDY 

With  respect  to  most  oceanic  basins,  in  the  Mediterranean  Sea  we  encounter  some 
unique  dynamic  and  physical  characteristics.  Such  features  as  the  vertical  structure, 
general  circulation,  deep  water  formation  and  strait  flows  over  the  sills  are  unique  to  this 
inland  sea.  The  Mediterranean  Sea  is  an  extreme  case  of  thermodynamic  interaction 
between  atmosphere  and  the  sea  with  the  sea  displaying  strong  reciprocal  responses  to 
the  air-sea  fluxes  of  mass  and  energy.  Because  of  the  large  heat  and  radiation  fluxes  into 
the  sea  and  excessive  evaporation,  the  Mediterranean  Sea  is  characterized  by  high 
salinities  and  temperatures,  relative  to  oceanic  water  masses  elsewhere. 

The  main  purpose  of  this  study  is  to  understand  the  effects  of  strong  evaporation 
and  other  seasonally  varying  atmospheric  forcing  on  mixed  layer  dynamics  in  the 
Mediterranean  Sea.  To  do  this,  the  mixed  layer  model  of  Garwood  (1977)  was  em- 
ployed, using  climatological  forcing  and  analyses  and  forecasts  from  Fleet  Numerical 
Oceanographic  Center  (FNOC).  Model  outputs  are  shown  in  Chapter  IV  in  this  thesis. 
In  Part  (IB),  the  statement  of  problem  is  described.  In  the  Part  (1C),  the  literature 
is  reviewed  and  summarized,  in  terms  of  the  dynamic  and  physical  characteristics  and  the 
general  circulation  of  the  Mediterranean  Sea. 

In  the  second  chapter  the  equations  on  which  the  OPBL  model  was  based  are  given. 
It  is  shown  how  from  those  equations  a  numerical  solution  is  obtained  to  solve  the  initial 
value  problem.  The  third  chapter  consists  of  observations  and  descriptive  discussion.  In 
the  fourth  chapter,  the  model  simulations  and  results  are  shown.  In  the  last  chapter, 
conclusions  and  recommendations  are  presented. 

B.  STATEMENT  OF  PROBLEM 

The  upper  layer  of  the  sea  responds  to  atmospheric  forcing  on  seasonal,  synoptic 
and  diurnal  time  scales  (Garwood,  1977).  In  this  study,  the  Naval  Postgraduate  School 
Oceanic  Planetary  Boundary  Layer  Model  (OPBL)  was  used  for  a  variety  of  cases  rep- 
resentative of  atmospheric  forcing  in  the  Mediterranean  Sea.  Such  cases  include  selective 
forcing  of  climatological  values  for  wind  stress,  net  evaporation  minus  precipitation 
(E-P)  rate,  and  heat  flux  values  to  demonstrate  individually  the  effect  of  seasonally 
varying  variables  on  the  Mediterranean  mixed  layer  dynamics.  The  model  was  applied 


at  three  stations  which  are  purposely  choosen  because  of  their  significant  differences  in 
atmospheric  forcing,  excess  evaporation  rate  and  general  circulation  with  respect  to  each 
other. 

In  the  Mediterranean  Sea,  convective  processes  play  a  very  important  role  in  the 
formation  of  Levantine  Intermediate  Water  and  deep  water.  In  order  to  understand  the 
water  mass  formation  process  within  the  Mediterranean  Sea,  we  need  to  know  the  nec- 
essary driving  mechanisms  including  wind  forcing,  heat  and  salt  fluxes,  thermohaline  ef- 
fects, vorticity  and  the  effects  of  the  general  circulation.  Due  to  the  dry,  cold  wind 
experienced  during  the  wintertime,  the  predominant  mixing  may  be  attiributed  to  con- 
vection cells  reaching  the  deeper  layers.  The  relative  roles  of  wind  stress  and  free  con- 
vection for  the  mixing  need  to  be  evaluated  quantitatively.  An  excess  of  evaporation 
over  precipitation  for  the  Mediterranean  Sea  requires  that  there  be  a  net  inflow  through 
the  Strait  of  Gibraltar  to  achieve  mass  balance  (Garrett,  1983). 

C.      LITERATURE  REVIEW 

1 .     General  Characteristics  of  The  Mediterranean  Sea 

The  vertical  structure  of  the  Mediterranean  Sea  consists  of  three  water  types: 

(i)  the  surface  Atlantic  water, 

(ii)  Levantine  Intermediate  Water,  and 

(iii)  the  Mediterranean  deep  water. 

Horizontally,  the  Mediterranean  has  two  geographical  regions:  the  western  and 
eastern  basins  which  are  divided  by  a  sill  of  400  m  depth  extending  from  Sicily  to  the 
North  African  coast.  The  maximum  depths  are  about  3400  m  in  the  western  basin  and 
4200  m  in  the  eastern.  The  Levantine  Intermediate  Water  is  formed  in  the  eastern  basin 
during  the  winter  off  the  south  coast  of  Turkey  with  a  temperature  of  15  °  C  and  a 
salinity  of  39.1  g/kg.  It  flows  westward  between  200  and  600  m  depth  along  the  North 
African  coast  and  out  below  the  surface  through  the  Strait  of  Gibraltar  into  the  Atlantic. 
In  the  Atlantic  this  intermediate  water  has  a  temperature  and  salinity  that  have  been 
modified  to  13  °  C  and  37.3  g/kg,  respectively,  by  mixing  with  overlying  Atlantic  water. 
The  outflow  of  Intermediate  water  below  the  surface  is  replaced  by  a  surface 
inflow  into  the  Mediterranean.  Even  though  the  western  Mediterranean  deep  water  ex- 
tends from  100  m  below  the  sill  of  Gibraltar  to  the  bottom,  it  is  surprising  to  find  that 
both  western  Mediterranean  deep  water  and  Levantine  Intermediate  Water  in  the 
Atlantic.  What  can  be  the  reasons  for  that?  Bryden  and  Stommel  (1982)  mention  the 
deep  westward  flow  along  the  Moroccan  continental  slope  with  an  anticyclonic  gyre 


deep  westward  flow  along  the  Moroccan  continental  slope  with  an  anticyclonic  gyre 
which  provides  the  high  velocities  necessary  to  raise  the  deep  water  up  toward  the  Strait 
of  Gibraltar.  Also,  high  velocities  in  the  Strait  could  provide  an  energy  source  to  raise 
water  from  great  depth  up  and  over  the  sill.  These  high  velocities  are  evident  in  the 
anticyclonic  gyre  in  the  Alboran  Sea  (Stommel,  1973). 

In  the  Mediterranean  Sea,  evaporation  exceeds  precipitation  plus  river  runoff. 
Therefore  there  is  a  net  loss  of  fresh  water  which  is  made  up  by  inflow  of  lesser  salinity 
Atlantic  water.  For  this  Atlantic  inflow  process,  the  Mediterranean  Sea  water  has  a 
residence  time  of  about  100  years. 

2.     Surface  Forcing  in  Mediterranean 

The  wind  stress  at  the  sea  surface  generates  turbulent  kinetic  energy  which  does 
work  against  the  buoyancy  forces,  entraining  the  underlying  denser  water  into  the  mixed 
layer  and  thus  increasing  the  potential  energy  of  the  water  column.  The  discontinuity  in 
temperature  and  density  that  occurs  at  the  bottom  of  the  mixed  layer  is  caused  by  the 
entrainment  process,  the  rate  of  which  is  inversely  proportional  to  the  magnitude  of  the 
density  jump. 

Changes  in  the  atmospheric  pressure  over  the  Mediterranean  are  the  main 
driving  mechanism  for  a  significant  part  of  low  frequency  sea  level  changes  (Garrett, 
1983).  The  wind-induced  set-up  outside  the  Strait  of  Gibraltar  is  less  important  than  is 
atmospheric  pressure  in  inducing  flows  through  the  Strait,  except  possibly  for  the  high 
frequency  response.  The  set-up  on  the  shelf  arises  from  a  northward  wind,  which  lags 
local  atmospheric  pressure  by  90  °.  So,  the  low  frequency  sea  level  changes,  and  Strait 
flows  are  found  to  be  correlated  with  large  scale  meteorological  forcing. 

Bryden  and  Stommel  (1984)  suggest  that  the  Mistral,  which  is  a  very  strong, 
cold,  dry  continental  wind  from  the  northwest  does  not  effect  the  deep  water,  unless  a 
very  strong  buoyancy  loss  associated  with  wintertime  convection  exceeds  a  certain  crit- 
ical value. 

According  to  the  POEM  (Physical  Oceanography  of  the  Eastern  Mediterranean) 
group's  report  (Eos  April  5,  1988),  the  eastern  Mediterranean  constitutes  a  general  cir- 
culating basin  in  itself,  including  important  wind  and  thermally  driven  components.  In 
this  basin,  the  convection  processes  produce  the  Levantine  Intermediate  Water  (LIW) 
in  the  return  flow  toward  Gibraltar.  This  flow  forms  the  main  body  of  the  North  Atlantic 
Upper  Deep  Water  at  intermediate  depth  in  the  Atlantic.  The  wind,  thermohaline  fluxes, 


It  is  suggested  that  the  Mistral  can  cause  the  convective  mixing  cells  capable  reaching 
the  intermediate  water  (LIW)  depth  as  well  as  deep  layers  in  the  southern  Adriatic.  They 
conclude  that  the  wind-driven  and  the  thermohaline  components  of  the  circulation  are 
predominately  driven  by  the  air-sea  latent  and  sensible  heat  fluxes  in  the  eastern 
Mediterranean  Sea  where  the  general  circulation  is  strongly  seasonal. 

In  terms  of  overall  mass  balance  it  is  expected  that  transport  into  the 
Mediterranean  should  exceed  transport  out  of  the  Mediterranean,  and  the  difference  is 
balanced  by  the  evaporation  process  over  the  entire  Mediterranean  Sea. 

3.     Deep  water  flow  and  formation  in  the  Mediterranean  Sea 

In  1970,  the  Medoc  Group  proposed  a  first  general  description  of 
Mediterranean  deep  water  formation,  divided  into  three  consecutive  phases: 

(i)  a  preconditioning  phase  which  is  the  erosion  of  the  upper  layer  vertical 

stratification; 

(ii)  strong  vertical  mixing  by  meteorological  forcing  (usually  caused  by  winds); 

and 

(iii)  sinking  and  spreading  due  to  a  baroclinic  instability,  producing  small  eddies. 
In  studying  the  formation  of  deep  water  in  the  northwestern  Mediterranean  Sea 
in  winter,  Gascard  (1978)  finds  the  frequent  cyclonic  and  anticyclonic  eddies  are  associ- 
ated with  processes  of  deep  water  formation.  He  suggests  that  a  baroclinic  instability  is 
responsible  for  the  slow  current  and  the  aperiodic  vertical  motions,  and  it  is  related  to 
the  formation  of  meanders,  eddies,  and  fronts  in  the  horizontal  plane.  Also,  the  north- 
west Mediterranean  basin  has  been  characterized  during  all  seasons  by  a  large  cyclonic, 
thermohaline  circulation  of  surface  water  masses  coming  from  the  Atlantic  and  subsur- 
face water  masses  originating  in  the  eastern  Mediterranean  (LIW)  above  a  thick  deep 
water  layer  (1500  to  2000  m).  The  Levantine  Intermediate  Water  has  been  formed  south 
of  Turkey  around  the  Islands  of  Rhodes  and  Cyprus  with  a  maximum  temperature  and 
salinity  value  at  a  depth  of  about  150-300  m.  During  the  winter  the  MEDOC  area 
(northwestern  Mediterranean  Sea)  is  subjected  to  the  effect  of  strong,  cold,  dry,  conti- 
nental winds  from  northwest  to  north  (Mistral  and  Tramontane).  High  evaporation  and 
strong  cooling  act  to  reduce  the  stability  of  the  central  part  of  the  surface  layer.  There 
is  a  hydraulic  control  on  the  flows  through  the  straits  in  combination  with  the  evapo- 
ration over  each  of  the  basins,  with  the  western  Mediterranean  surface  water  being  1.7 
g/kg  saltier  than  the  Atlantic  surface  water,  and  the  eastern  Mediterranean  0.6  g/kg 


g  kg  saltier  than  the  Aiantic  surface  water,  and  the  eastern  Mediterranean  0.6  g'kg 
saltier  than  the  western  Mediterranean.  The  deep  water  in  the  eastern  Mediterranean 
appears  to  be  formed  in  the  Adriatic  (Pollak,  1951),  which  is  separated  from  the  eastern 
basin  by  the  Strait  of  Otronto.  The  MEDOC  group  remark  that  the  deep  water  forma- 
tion in  the  eastern  basin  is  not  a  late-winter  conversion  of  intermediate  water  into  deep 
water,  because  the  Adriatic  is  a  separate  basin  with  little  intermediate  water  in  its 
northernmost  part. 

According  to  Millot  (1986)  the  fractional  amount  of  water  of  Atlantic  origin 
involved  in  the  deep  water  formation  processes  occurring  south  of  the  French  coast  is 
roughly  one-fourth  (the  other  three-fourths  being  Levantine  Intermediate  Water)  of  the 
volume  of  water  exiting  from  the  Mediterranean  Sea  through  the  Strait  of  Gibraltar.  So, 
the  flow  of  LIW,  although  disturbed  by  various  mesoscale  phenomena,  follows  a  steady 
cyclonic  path  along  the  continental  slope  off  Italy,  France  and  Spain,  and  it  appears  to 
be  mainly  forced  by  buoyancy.  Deep  water  formation  occurs  during  winter  in  the  Gulf 
of  Lions  (Gascard,  1978).  The  reasons  for  the  circulation  inside  the  western 
Mediteranean  Sea  are  not  well  understood  at  present.  Lacombe  et  a/.(1980)  show  that 
the  density  stratification  in  the  Mediterranean  Sea  is  controlled  during  the  May-October 
warm  season  by  the  thermal  structure.  During  the  cold  season  it  is  controlled  mainly  by 
the  salinity  structure.  In  summer  time  the  mean  vertical  motion  is  very  weak.  So  the  lo- 
cal deep  water  formation  should  occur  during  downward  convective  events  of  short  du- 
ration. In  the  presence  of  large  thermal  losses  in  winter,  due  to  high  winds  (Mistral  and 
Tramontane)  and  dry  cold  air  blowing  over  a  warm  sea,  the  density  reached  in  the  cen- 
tral area  is  greater  than  elsewhere,  and  convective  motions  are  likely  to  start  there  first. 
This  is  called  the  preconditioning  phase.  During  a  strong  mixing  phase,  high  winds  last- 
ing for  ten  days  and  resulting  in  an  evaporation  rate  of  about  2  cm/day  are  able  to 
homogenize  the  water  down  to  2000  m. 

4.     Mixing  and  Vertical  Convection  Processes  in  the  Mediterranean 

As  in  most  temperate  seas,  the  depth  of  the  surface  mixed  layer  is  typically  on 
the  order  of  ten  to  one  hundred  meters,  and  the  horizontal  scale  size  is  that  of  the 
internal  Rossby  radius,  typically  12  to  15  km  in  the  Mediterranean  Sea. 

Deepening  of  the  mixed  layer  is  caused  by  entrainment  of  the  more  dense 
underlying  water  into  the  turbulent  region  above.  This  process  leads  to  a  potential  en- 
ergy increase  and  cannot  take  place  without  the  turbulent  kinetic  energy  generated  by 


the  wind  stress  and  buoyancy  loss  at  the  surface  of  the  mixed  layer.  Progressive  cooling 
and  evaporation  which  occur  permanently  in  winter  during  periods  of  calm  act  on  the 
quasi-homogeneous  mixing  zone  and  may  induce  a  slow  downward  mean  motion  of  new 
homogeneous  and  dense  waters.  This  is  the  first  process  of  deep  water  formation,  it  is 
called  the  advective  process. 

Also,  strong  and  sudden  surface  cooling  and  evaporation  produced  by  the 
Mistral  or  Tramontane  acting  on  such  a  homogeneous  and  unstable  structure  may  trig- 
ger intense  deep  conveciion  within  the  interior  of  the  fluid.  This  is  called  the  connective 
process,  and  it  causes  the  formation  of  new  deep  waters.  A  principal  goal  of  this  thesis 
research  is  to  quantify  the  relative  importance  of  wind-driven  forced  convection  and 
buoyancy-driven  free  convection  during  this  convective  stage. 

According  to  the  Pollak  (1951)  the  vertical  convection  currents  can  cause  a 
complete  overturning  of  the  entire  water  column  and  can  thus  be  an  important  factor  in 
the  formation  of  subsurface  water  characteristics  which  are  most  likely  to  occur  during 
periods  of  minimum  surface  temperatures.  Pollak  also  suggested  that  the  Adriatic  Sea, 
the  only  other  possible  source  of  deep  eastern  Mediterranean  water,  has  a  complicated 
structure.  It  is  highly  affected  by  the  short-term  weather  fluctuations  due  to  its  restricted 
size  and  mean  depth.  Under  those  conditions  the  mixed  water  sinks  and  then  flows 
sporadically  along  the  bottom  into  the  Mediterranean  basin,  where  it  circulates 
counterclockwise. 


II.      THEORY 

A.  OCEANIC  PLANETARY  BOUNDARY  LAYER  (OPBL)  MODEL 
1.     General  characteristics  and  dynamics  of  the  OPBL 

The  OPBL  model  is  a  one-dimensional,  second  order  turbulence  closure,  verti- 
cally integrated  (bulk)  model  of  the  upper  ocean  surface  turbulent  boundary  layer  de- 
veloped by  Garwood  (1977).  It  uses  the  continuity  equation  for  an  incompressible  fluid, 
the  first  law  of  thermodynamics  (heat  equation),  the  conservation  of  salt  equation,  an 
analytical  equation  of  state,  the  Navier-Stokes  equation  of  motion  with  the  geostrophic 
component  eliminated,  and  a  two-component  vertically  integrated  turbulent  kinetic  en- 
ergy (TKE)  budget. 

The  oceanic  boundary  layer  or  mixed  layer  is  the  fully  turbulent  region  that  is 
bounded  above  by  the  cir-ocean  interface,  and  it  is  where  the  temperature  and  salinity 
are  usually  observed  to  be  fairly  well  mixed.  Below,  the  mixed  layer  is  assumed  to  be 
bounded  by  a  dynamically  stable  thermocline.  It  is  assumed  that  the  turbulence  of  the 
overlying  mixed  layer  provides  the  energy  needed  to  destabilize  and  erode  the  underlying 
stable  thermocline  in  accordance  with  Garwood  (1977).  The  wind  and  upward  surface 
buoyancy  flux  through  the  surface  are  the  sources  of  mechanical  energy  for  the  gener- 
ation of  this  turbulence. 

The  model  differs  from  earlier  models:  First,  the  amount  of  wind-generated 
TKE  to  be  used  in  mixing  is  a  function  of  the  ratio  of  the  mixed  layer  depth  to  the 
Obukhov  mixing  length  L.  Second,  viscous  dissipation  is  dependent  on  a  local  Rossby 
number.  Finally,  separate  vertical  and  horizontal  equations  for  TKE  are  used,  permit- 
ting a  more  consistent  interpretation  of  both  entraining  and  retreating  mixed  layers. 

This  model  shall  be  used  to  help  understand  the  relative  roles  of  wind  stirring 
and  surface  buoyancy  fluxes  (attributable  to  both  temperature  and  salinity  fluxes)  in  the 
deep  water  formation  for  the  Mediterranean  Sea.  In  the  following  sections  it  is  shown 
how  to  include  salinity  and  the  effects  of  strong  evaporation  in  the  theory  for  application 
to  the  Mediterranean  Sea. 

B.  THE  OPBL  MODEL  EQUATIONS 

1.      Temperature  and  Salinity  Equations 

Application  of  the  first  law  of  thermodynamics  to  mixed  layer  entrainment 
gives, 


3T=-ATShjh.  (2.1) 

Where  ST  is  the  mixed  layer  temperature  reduction  due  to  entrainment,  AT  is  the 
discontinuity  between  the  mixed  layer  temperature  and  the  temperature  immediately 
below  the  mixed  layer.    The  mixed  layer  depth  h  is  the  distance  over  which  turbulent 
energy  must  be  transported  by  the  vertical  component  of  turbulent  velocity. 
The  mixed  layer  temperature  equation  is  therefore, 

_d]L__ALllL  m\ 

dt  "      h    dt  ■  Kl-l) 

The  time  derivative  of  h, 

^-W.-»HB  (2-3) 

is  the  entrainment  velocity  which  depends  on  the  vertical  velocity, 

-§7=^  (2.4) 

where  W(_h)  =  0. 

Then  the  entrainment  heat  flux  is 

-ATWe=TW{_hy  (2.5) 

For  mixed  layer  salinity,  conservation  of  salt  mass  in  the  presence  of 
entrainment  yields  a  set  of  equations  that  are  mathematically  analogous  to  equation 
(2.1),  equation  (2.2)  and  equation  (2.5): 

SS  =  -  ASdhjh.  (2.6) 

dS      _  AS    8h  (~  -x 

dt  ~      h    dt  •  {    } 

AS^  =  5W(_ft).  (2.8) 

Here  55  is  the  change  in  mixed  layer  salinity  due  to  an  entrainment  increase 
of  h  by  the  depth  increment  6h,  with  AS  being  the  salinity  drop  across  the  bottom  of  the 


mixed  layer.  In  equation  (2.8),    S'w\_h)  is  the  entrapment  salinity  flux,  and  IV,  is  the 
same  entrainment  velocity  invoked  in  equation  (2.3)  and  equation  (2.4). 

Including  the  surface  heat  flux,  equation  (2.2)  is  augmented  to  become, 

eT       (Q0lPCp-ATH't) 

ct  h  *  "  ' 

The  mixed  layer  salinity  budget  equation  (2.7)  is  similarly  augmented  to  include 
the  surface  salinity  flux,  S(E-P),  where  E-P  is  the  net  evaporation  minus  precipitation 
rate: 

es       S(E  -  P)  -  WAS 

-^-  =  — -, £=-.  (2.10) 

ct  h 

2.     Buoyancy  and  momentum  equations 

The  buoyancy  equation  is  generated  by  the  combination  of  the  heat  and  salt 
equations  with  an  equation  of  state, 

p  =  p0[l-a(Ar)  +  /?(AS)].  (2.11) 


where  y.  and  /?  are  the  expansion  coefficients  for  heat  and  salt. 
Buoyancy  is  defined  as, 

b  =  S      n,      .  (2-12) 


where  b  =  B  +  b',  decomposing  into  mean  and  turbulent  parts. 

The  mean  buoyancy  equation  is  therefore  derived  from  the  heat  and  salt 
equations  (2.9)  and  (2.10),  using  the  equation  of  state  (2.11)  and  the  definition  for 
buoyancy  (2.12)  : 

ct  cz  p0Cp 

By  assuming  the  Boussinesq  approximation,  no  viscous  terms,  incompressible 
fluid  and  subtracting  the  geostrophic  equations  from  the  total  momentum  equations, 
we  obtain  the  mean  momentum  equations: 


lLL=jy-Il^L  (2.14) 

&L-fU-&jL.  (2.15) 

at  oz 

Integration  of  (2.13),  (2.14),  and  (2.15)  over  the  entrainment  zone  gives  us  the 
jump  conditions  for  the  entrainment  turbulent  fluxes  of  buoyancy  and  momentum  re- 
spectively: 

-bW{_h)  =  ABJt  (2.16) 


■wy-fr  (2.i7) 


■-AV~-  (2.18) 


where 


dt 

b'  =  agT-fJgS' 
AB  =  agAT  -  pgAS 
AT=f-T(_h_6) 
AS  =  S-S{_h_6). 


Here  AT  and  AS  are  temperature  and  salinity  discontinuity  at  the  mixed  layer  base. 

3.     Mean  Turbulent  Kinetic  Energy  Equation 

By   subtracting   the   scalar   product   of  (u,v,w)  with   the   mean   momentum 
equations  from  the  mean  kinetic  energy  equation  we  get  the  mean  TKE  equation: 


|^  =  _[wV(^)  +  vV(^L)]  +  ,w_^K(|.  +  ^)]_^0.(2.19) 

(I)  (ID  (HI)  (IV) 

where 


E=u'z  +  v'z  +  w'\  (2.20) 

The  mixed  layer  turbulance  is  assumed  to  be  in  a  quasi-steady  state  because  the 
time  rate  of  change  TKE  is  usually  smaller  than  the  other  terms  and  will  be  neglected. 
The  term  (I)  represents  the  rate  of  mechanical  production  and  is  the  dominant  source 
of  TKE  for  wind  driven  regimes.  The  term  (II)  represents  the  buoyancy  flux  and  can 
be  either  a  source  or  a  sink.  For  strong  evaporation  and  net  surface  heat  loss  b'W  will 
be  large  and  positive  and  may  become  even  more  important  than  term  (I)  in  producing 
mixed  layer  turbulence.  The  term  (III)  is  the  divergence  of  the  turbulent  flux  of  TKE 
and  is  responsible  for  transporting  turbulence  generated  near  the  surface  to  the  bottom 
of  the  mixed  layer  where  it  can  cause  entrainment  and  mixed  layer  deepening.  The  term 
(IV)  represents  viscous  dissipation,  the  major  loss  of  TKE. 

The  TKE  equation  in  the  entrainment  zone,  z  =  -h  is  reduced  to  : 

jfi^-U  =  (  -*g*T+0gAS)Wt  +  Jtf*  -^p-  .  (2.21) 

II  III 

In  the  entrainment  zone,  E  generated  higher  in  the  layer  is  transported  (by  term 
III)  into  contact  with  the  underlying  water.  When  the  TKE  mixes  this  underlying  water 
with  the  mixed  layer  water,  at  a  rate  lVe  ,  the  TKE  is  damped  by  the  entrainment 
buoyancy  flux  (term  II).  In  the  entrainment  zone  (terms  I  and  IV)  are  neglected,  with 
the  dominant  balance  assumed  to  be  transport  and  buoyant  damping. 

4.     Bulk  Turbulence  Equations 

The  resultant  OPBL  model  equations  may  be  summarized  as  follows: 
Mixed  layer  depth: 

^=We-W{_hy  (2.22) 

Entrainment  rate,  IV„  is  derived  from  equation  (2.21): 


Vir»<£> 

e      hagAT+<E> 

Equation  (2.19)  is  integrated  vertically  to  give  total  TKE  for  the  mixed  layer: 


-(h<E>)  =  G0  +  Gh-B0-Bh-D.  (2.24) 


The  Vertical  TKE  component  becomes  : 

-L{h<W2>)  =  -B0-Bh  +  R-^-.  (2.25) 

Mixed  layer  currents  are,  equation  (2.26)  and  equation  (2.27): 

JL{h<u>)=fh<v>  +|,  (2.26) 


A-{h<v>)=fh<u>  +-f.  (2.27) 

Without  salinity,  the  buoyant  damping 'production  terms  at  the  surface  and  at 
the  bottom  are: 

B0  =  aghTW(0)  =  -aighQoKpCp).  (2.28) 


Bh  =  aghATWe.  (2.29) 

With  salinity  effect,  downward  surface  buoyancy  flux,  AB  and  Bh  become  re- 
spectively: 

*o  =  l«glQ0l(pCp)-]  -  pgS(E  -  P)-]h,  (2.30) 

AB  =  agAT-jlgAS,  (2.31) 

Bh  =  {*gAT-figAS)We.  (2.32) 

The  shear  production  at  the  surface  and  at  the  bottom  are  : 

G0=12tf3  =  12[(^)2  +  (^)2]3/4  (2.33) 

Gh  =  [<U>2+  <v>2We-  (2-34) 


Dissipation  is  modeled  as, 

Z)=2<£>3/2.  (2.35) 

Redistribution,  R  in  (2.25)  is, 

R  =  2{<E>  -3<w!2>)yfE  (2.36) 

The  following  are  the  constants  and  the  parameters  which  are  needed  in  the 


model: 


a.  Model  Constants 

pc    Density  of  air 

p0    Density  of  sea  water 

a     Thermal  expansion  coefficient 

/?     Haline  expansion  coefficient 

g    Acceleration  of  gravity 

Cp    Specific  heat  at  constant  pressure 

CD     Drag  coefficient 

b.  Model  Parameters 

f    Coriolis  parameter 

t     Wind  Stress 

QQ     Incoming  solar  radiation 

Qb    Longwave  back  radiation 

Q,    Latent  heat  flux 

Qh    Sensible  heat  flux 

E    Evaporation  rate 

P    Precipitation  rate 

5.      Approximation  For  Entrainment  Equation 

a.     Preliminary  Assessment  of  Effects  of  Evaporation 

The  Garwood  (1977)  model  system  equations  (2.22)  to  (2.36)  may  be  ap- 
proximated by  the  Krauss-Tumer  (1967)  model,  which  neglects  dissipation  (D)  and 
entrainment  shear  production  (Gh): 


We  =  (GQ-B0)IAB  (2.37) 

b.      Steady-State  Non- Entraining  Mixed  Layer  (Spring  and  Summer) 

If   B0    is    positive    and    equals    G0,    then    there    is    no    entrainment, 
W,  =  dhjdt  =  0,  and  B0  =  G0  ,  or 

[ag(  ■&.  C,)  -  &S(£  -  P)]A  =  G0.  (2.38) 

Solving  for  the  steady-state  diagnostic  mixed  layer  depth  in  equation  (2.38) 
gives: 


l«gQolpCp-PgS(E-P)] 


(2.39) 


It  is  readily  seen  that  a  net  (E—P)>0  will  reduce  the  denominator,  in- 
creasing h.  In  addition  to  an  obvious  increase  in  mixed  layer  salinity  due  to  evaporation, 
the  mixed  layer  temperature  will  be  indirectly  changed  by  the  increase  in  h.  So  this  latter 
effect,  equation  (2.8)  is  combined  with  equation  (2.39)  with  We  =  0, 

n Q^c^ 

ct        IGolagQolpCp-PgSiE-P)! 

In  the  numerator  of  equation  (2.40)  evaporation  decreases  Q0  directly  due 
to  the  loss  of  latent  heat.  But,  the  accompaning  increase  in  mixed  layer  depth  will  further 
reduce  mixed  layer  temperature  by  increasing  the  denominator  of  the  right-hand  side  of 
equation  (2.40).  Thus  it  has  been  shown  that  we  should  expect  the  surface  salinity  flux 
due  to  E  to  indirectly  influence  the  evolution  of  the  temperature  profile,  as  well  as  the 
evolution  of  the  salinity  profile. 

c.      Deepening  Mixed  Layer  (Fall  and  Winter) 

During  the  fall  and  winter,  the  mixed  layer  deepens  and  W,  #  0.  Therefore 

the  complete  model  equations  must  be  solved.    Nevertheless,  the  effects  of  W>§  may 

be  seen  in  the  heat  and  salt  equations  (2.6)  and  (2.10).  The  general  tendency  for  change 

in  mixed  layer  temperature  due  to  a  strong  evaporation  will  be  due  to  two  contributions: 

(i)  <20  will  be  made  more  negative  due  to  the  loss  of  latent  heat, 


(ii)  h  and  U'e  will  be  increased  due  to  the  increased  production  of  TKE  by 
the  extra  surface  buoyancy  flux.  In  the  Mediterranean  AT  will  almost  always  be  posi- 
tive, so  enhanced  entrainment  will  always  reduce  the  mixed  layer  temperature. 

In  the  mixed  layer  salinity  budget,  equation  (2.10),  the  net  effect  of  strongly 
positive  (E-P)  is  to  increase  mixed  layer  salinity.  Whereas,  equation  (2.10)  shows  that  the 
increase  is  not  as  great  as  it  will  be  for  a  constant-depth  mixed  layer  (as  in  the  summer). 
The  —  W£S  term  of  equation  (2.10)  will  usually  reduce  mixed  layer  salinity  because  AS 
will  generally  be  positive  during  the  fall  and  winter.  However,  it  may  be  negative  in  some 
locations,  especially  if  there  is  strong  horizontal  advection  with  higher  salinity,  water 
interleaving  below  the  mixed  layer.  Also,  (£-  P)>  0  increases  h,  reducing  the  magnitude 
of  dS/dt. 


III.      OBSERVATIONS  AND  DESCRIPTIVE  DISCUSSION 


A.      OBSERVED  SURFACE  HEAT  FLUX  CHARACTERISTICS 

In  our  problem  daily  heat  and  radiation  fluxes  are  provided  from  the  climatology 
of  Hsiung  (1983,  1985)  for  the  three  locations  of  interest  in  the  Mediterranean  Sea. 
The  net  surface  heat  flux  can  be  written  as, 

Qo  =  Qs-Qe-QH~Qb  (4-1) 

where 

Qs  Incoming  Solar  Radiation, 

Oe  Upward  Latent  Heat  Flux 

Qh  Upward  Sensible  Heat  Flux 

Qb  Outgoing  Longwave  Radiation 

The  latent  heat  flux  proportional  to  E, 

Qe  =  pLE  (4.2) 

where  the  evaporation  rate  E  in  cm/ day  is  also  available  from  climatological  data. 

The  fluxes  of  radiation,  latent  and  sensible  heat  at  the  sea  surface  determines  shal- 
lowing or  deepening  of  the  mixed  layer  by  entrainment.  According  to  the  equations 
(2.33),  (2.36)  and  (2.39),  if  there  is  a  positive  buoyancy  flux  {Q0)<0  and  (E>  P)  then  the 
mixed  layer  must  deepen.  This  will  be  the  case  at  night  or  whenever  longwave  radiative 
cooling  plus  the  upward  turbulent  fluxes  of  heat  exceed  solar  radiation.  The  negative 
downward  buoyancy  flux  will  result  in  the  shallowing  of  the  MLD  provided  wind  stress 
is  not  large. 

In  the  Mediterranean,  with  large  (E-P),  it  may  be  possible  to  have  a  net  downward 
heat  flux,  Q0>0,  and  still  have  a  net  upward  surface  buoyancy  flux.  From  equation 
(2.30),  B0<0  whenever 

l*g-f^--feS(E-P)l<0 


Solving  the  inequality  for  E-P,  this  situation  will  occur  if, 


(£-/>)>[— ^p-].  (4.3) 

This  effect  may  be  particularly  important  during  the  late  winter  or  early  spring 
when  O0  may  be  positive,  but  small.  If  this  were  to  happen,  the  surface  buoyancy  flux 
would  be  positive,  producing  added  TKE  to  possibly  deepen  the  mixed  layer  well  beyond 
the  maximum  depth  possible  for  a  temperature-dominant  ocean.  This  process  could  also 
be  an  important  mechanism  to  mix  heat  downward,  overcoming  the  normal  stratifying 
tendency  of  the  surface  heat  flux. 

The  distribution  of  the  surface  heat  flux  does  not  have  significant  differences  over 
Mediterranean  from  place  to  place.  However,  the  large  evaporation  rate  over  the 
Mediterranean  is  very  unusual  and  it  causes  the  strong  turbulent  mixing  in  the  surface 
boundary  layer.  Due  to  the  similar  surface  heat  flux  values  over  the  Mediterranean  Sea 
we  do  not  expect  significantly  different  results  in  terms  of  the  mixed  layer  dynamics  be- 
tween three  stations  which  are  chosen  in  the  Mediterranean  Sea. 

From  Fig.  1  on  page  18,  Fig.  2  on  page  19,  Fig.  3  on  page  20,  and  Fig.  4  on  page 
21,  we  see  the  latent  heat  flux  and  incoming  solar  radiation  values  at  two  stations  in  the 
central  and  western  Mediterranean,  respectively. 

At  the  eastern  Mediterranean  Station  [<f>  =  35°AT,  X  =  30°£),  the  latent  heat  flux 
decreases  from  the  initial  value  which  is  165  W/m2  to  43  W/m2  in  winter.  After  May,  it 
increases  until  end  of  the  year.  After  the  middle  of  October,  the  latent  heat  flux  increases 
sharply,  probably  due  to  the  large  air-sea  temperature  difference.  Depending  on  the  la- 
tent heat  flux,  evaporation  contributes  to  the  buoyant  production  of  TKE  and  causes 
vertical  mixing.  Therefore,  at  this  station  the  TKE  is  enhanced  by  evaporation  the  most 
during  November. 

Incoming  solar  radiation  is  almost  a  sinusoidal  curve  over  the  year.  See  Fig.  5  on 
page  22.  It  increases  smoothly  from  the  January  value  81  W/m2  until  July,  then  decreases 
smoothly  until  end  of  the  year. 

Fig.  7  on  page  24  shows  the  net  heat  flux.  From  January  to  July  it  increases  from 
-183  W  m2  to  211  W/m2,  and  it  then  decreases  to  -200  W/m2  in  December. 

At  the  central  Mediterranean  station  ((f)  =  35° N,  X  =  20°£)  ,  Fig.  2  on  page  19  the 
latent  heat  flux  differs  from  that  at  the  eastern  station.  At  this  station,  latent  heat  flux 
peaks  in  February,  then  decreases  until  the  middle  of  the  June.  It  then  increases  through 
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Fig.   1.      The  Latent  Heat  Flux  At  40  °  N  05  °  E 


October,  with  an  unusual  plateau  in  August  and  September.  After  October  it  decreases 
to  144  W/m2. 

Incoming  solar  radiation  in  central  part  of  Mediterranean  differs  from  in  the  eastern 
part  only  slightly.  Also  the  net  downward  heat  flux  is  similar  to  that  in  the  eastern 
Mediterranean. 

In  the  western  Mediterranean,  Fig.  1  shows  the  latent  heat  flux.  From  January  to 
May,  latent  heat  flux  and,  therefore,  evaporation  rate  gets  smaller.  It  then  increases 
slightly  in  June.  Between  June  and  July  it  decreases  slightly.  From  July  until  November, 
latent  heat  flux  increases  sharply,  then  decreases  until  end  of  the  year.  The  maximum 
latent  heat  flux  is  in  November.  Minima  are  in  May  and  July. 
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Fig.  2.      The  Latent  Heat  Flux  At  35  °  N  20  °  E 


Incoming  solar  radiation  is  similar  to  the  other  stations.  See  Fig.  3  on  page  20, 
Fig.  4  on  page  21  and  Fig.  5  on  page  22. 

B.      THE  OBSERVED  WIND  STRESS  DISTRIBUTION 

The  wind  stress  values  are  computed  from  climatological  wind  speed  data.  The  drag 
coefficient,  CD  depends  upon  the  wind  speed  climatology  and  the  formulation  of  Bunker 
(1976). 

As  for  the  surface  heat  fluxes,  the  wind  velocity  values  don't  differ  significantly  for 
the  different  locations.  However,  the  western  station  has  somewhat  stronger  wind 
speeds  than  does  the  eastern  Mediterranean.  See  Fig.  8  on  page  25,  Fig.  9  on  page  26, 
and  Fig.   10  on  page  27. 
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Fig.  3.      The  Incoming  Solar  Radiation  At  40  °  N  05  °  E 


From  seasonal  variations  in  the  wind,  we  subdivide  the  year  into  three  periods.  In 
the  first  period,  (between  January  and  June  in  western  and  central  Mediterranean,  but 
from  January  to  May  in  eastern  Mediterranean),  wind  stress  generally  is  decreasing.  In 
the  second  period  which  is  in  summer,  wind  speeds  are  lowest.  This  period  lasts  usually 
until  September.  Then  in  the  third  portion  of  the  year,  wind  stress  increases  continually. 
So,  the  maximum  wind  stress  occurs  between  November  and  April,  and  the  minimum 
occurs  between  May  and  October. 

Fig.  8  on  page  25,  Fig.  9  on  page  26  and  Fig.  10  on  page  27  show  the  wind  velocity 
as  functions  of  time  for  the  selected  stations  in  the  Mediterranean  Sea. 


Fig.  4.      The  Incoming  Solar  Radiation  At  35  °  N  20  °  E 


C.      OBSERVED  TEMPERATURE  AND  SALINITY  CHARACTERISTICS 

The  eastern  Mediterranean  basin  is  isolated  from  the  western  basin  by  the  Strait 
of  Sicily.  The  Eastern  basin  experiences  some  important  wind  and  thermal  forcing.  Due 
to  convection  processes,  the  salty  Levantine  Intermediate  Water  is  produced  here.  This 
water  mass  (LIW)  flows  into  the  Atlantic  Ocean,  and  it  forms  the  main  body  of  the 
North  Atlantic  upper  deep  water.  (Malonette  et  ai,  1988).  Understanding  the  forma- 
tion of  this  water  mass  is  therefore  important  for  the  Atlantic,  as  well  as  the  local 
thermohaline  circulations.  The  eastern  Mediterranean  has  a  maximum  depth,  surface 
temperature  and  salinity.  4200  m.  26  °  C  and  39.1  g  kg.,  respectively.  Fig.  11  on  page 
28  shows  the  observed  "emperature  structure  in  the  eastern  Mediterranean  Sea.   As  can 


Fig.  5.      The  Incoming  Solar  Radiation  At  35  °  N  30  °  E 


be  seen  from  the  complete  picture,  the  oceanographic  year  can  be  considered  as  having 
a  cool  season  from  November  to  April  and  warm  season  from  May  to  October.  During 
the  warm  season,  mixed  layer  temperature  increases,  and  the  16.5  °  C  isotherm  rises  to 
170  m  depth  due  to  the  surface  heating  and  calm  weather  conditions.  Below  the  170  m 
depth  the  vertical  structure  is  most  stratified.  In  the  cool  season,  the  surface  layer  expe- 
riences strong  wind  effects  and  surface  cooling  due  to  atmospheric  winter  conditions  and 
excessive  evaporation  rate.  Therefore,  buoyancy  loss  causes  vertical  mixing  within  the 
surface  boundary  layer,  and  this  mixing  can  penetrate  to  a  depth  of  1200  m.  See 
Fig.  12  on  page  29. 


Fig.  6.      The  Latent  Heat  Flux  At  35  °  N  30  °  E 


At  the  central  Mediterranean  Sea  station,  ((f)  =  35°A",  k  =  20°£)  the  mixed  layer 
temperature  drops  due  to  winter  atmospheric  conditions.  The  MLD  deepens  below  100 
m.  In  summer,  due  to  calm  weather  conditions  and  excess  downward  heat  flux,  mixed 
layer  temperature  increases  and  is  accompanied  by  a  significant  reduction  in  MLD. 

The  most  significant  differences  between  the  central  and  the  eastern  Mediterranean 
temperature  structure  are  attributable  to  the  differences  in  surface  heat  flux  in  the  east- 
ern Mediterranean  relative  to  the  surface  heat  flux  in  the  central  regime.  There  are  higher 
temperatures  at  all  depths  in  the  eastern  regime  relative  to  the  central.  See  Fig.  12  on 
page  29  and  Fig.   13  on  page  30.  Another  significant  difference  is  that  the  wind  stress  in 
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Fig.  7.      The  Net  Heat  Flux  At  35  °  N  30  °  E 


the  central  regime  is  stronger  than  in  the  eastern  regime  in  any  month  of  the  year.  These 
two  atmospheric  forcing  factors  play  a  very  important  role  in  the  mixed  layer  dynamics. 

From  Fig.  14  on  page  31,  the  surface  salinity  value  is  highest  between  May  and 
October,  probably  due  to  relatively  warmer  and  calm  wind  conditions,  coupled  with 
higher  evaporation  rate. 

The  shallow  summer  time  mixed  layer  is  responsible  for  this  increase  in  surface 
salinity.  The  strong  evaporation  continues  in  spite  light  winds.  However,  the  light  winds 
(small  G0)  and  strong  insolation  (large  Q0)  cause  the  mixed  layer  to  be  shallow,  concen- 
trating the  salinity  in  a  small  volume.  This  is  predicted  by  equation  (2.24)  : 
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Fig.  8.      The  Wind  Velocity  At  40  °  N  05  °  E 
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Here  £  -  P  is  insufficiently  large  to  overcome  Q0  ,  and  the  inequahty  (4.3)  is  not 
yet  satisfied  : 


(E-P)<  <ag 


pPcps 


However,  as  fall  comes  on  and  QQ  is  reduced,  the  dominant  buyancy  flux  contrib- 
ution shifts  from  Q0  to  (E-P).    This  salty  but  warm  surface  water  cools  in  fall  and  be- 


Fig.  9.      The  Wind  Velocity  At  35  °  N  20  °  E 


comes  denser  than  the  water  below,  and  it  sinks.  Fig.  14  on  page  31  shows  an  example 
of  convected  surface  water  with  a  core  depth  of  210  m  and  160  m.  As  is  suggested  by 
Malanotte-Rizzoli  and  Robinson  (1987),  the  Mistral  and  Tramontone  winds  which  are 
strong,  cold  and  dry  can  start  the  mixing  process,  called  the  preconditioning  phase  .  This 
can  overturn  the  water  column  down  to  a  couple  of  hundred  meters  depth.  This  process 
is  equivalent  to  the  formation  process  of  Levantine  Intermediate  Water  in  the  eastern 
Mediterranean  which  was  also  suggested  by  Bryden  and  Stommel  (1983). 

Unlike  the  eastern  Mediterranean,  salinity  increases  with  depth  in  the  central 
Mediterranean.  See  Fig.  15  on  page  32  and  Fig.  16  on  page  33.  This  reduction  in  tem- 
perature and  increase  in  salinity  with  depth  makes  the  vertical  water  structure  very  stable 


Fig.   10.      The  Wind  Velocity  At  35  °  N  30  °  E 


in  this  pan  of  the  Mediterranean  Sea.  Due  to  such  a  stable  vertical  structure  in  the 
central  Mediterranean  Sea,  the  mixing  process  can  only  be  significant  during  the  strong 
wind  storms  in  winter  associated  with  the  Mistral  and  Tramontone. 


OBSERVED  TEMPERRTURE  P.T  35  N     3D  E 


Fig.   1 1.      The  Observed  Temperature  At  35  °  N  30  °  E 
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Fig.   12.      The  Observed  Temperature  At  35  °  N  30  °  E 
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Fig.   13.      The  Observed  Temperature  At  35  °  N  20  °  E 
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Fig.   14.      The  Observed  Salinity  At  35  °  N  30  °  E 
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Fig.   15.      The  Observed  Salinity  At  35  °  N  20  °  E 
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Fig.   16.      The  Observed  Salinity  At  35  °  N  20  °  E 


IV.      OPBL  MODEL  SIMULATIONS 

The  OPBL  Model  is  used  to  study  the  relative  importance  of  evaporation,  heat  flux, 
and  wind  stress  for  mixing  at  three  stations  in  the  Mediterranean  Sea.  The  monthly  heat 
flux  data  for  these  Mediterranean  locations  were  obtained  from  the  climatological  data 
computed  by  Hsiung  (1983,  1985),  and  the  wind  stress  data  were  provided  from  wind 
stress  climatology.  (U.S.  Navy  Marine  Climatic  Atlas).  The  specified  daily  boundary 
conditions  are  latent  heat  flux  (Qt)  ,  incoming  solar  radiation  ((?,),  net  heat  flux  (Q,)  , 
wind  stress  (t),  evaporation  (E),  and  precipitation  (P).  The  predicted  model  variables  are 
mixed  layer  depth  (h),  temperature  (T),  salinity  (S),  and  below-layer  profiles  T(z)  and 
S(z). 

The  model  is  run  at  the  three  locations,  using  three  cases: 

(i)  free  convection  (heat  and  salinity  fluxes)  only  ; 

(ii)  forced  convection  plus  heat  flux  only  ;  and 

(iii)  combined  free  and  forced  convection  (heat,  salinity  and  momentum  fluxes) 
A  discussion  of  model  results  follows. 

A.     CASE  1   (FREE  CONVECTION  ONLY,  NO  WIND  STIRRING) 

The  purpose  of  this  case  is  to  show  the  effect  of  excessive  evaporation  over  pre- 
cipitation without  wind  effects  on  the  mixed  layer  dynamics  at  all  three  stations.  One 
can  see  the  correlation  between  net  heat  flux  Fig.  17  on  page  35  and  mixed  layer  depth 
Fig.   18  on  page  36. 

For  the  case  where  the  evaporation  exceeds  precipitation  in  the  eastern 
Mediterranean,  creating  a  net  evaporation,  the  mixed  layer  temperature  changes  depend 
on  the  seasonally  varying  net  heat  flux.  In  January  mixed  layer  depth  deepens  from  in- 
itial depth  (95  m)  to  197  m.  After  January  mixed  layer  depth  is  constant  at  197  m  until 
the  second  half  of  April.  At  the  end  of  April  the  mixed  layer  shallows  rapidly  to  0  m,  and 
sudden  fluctuations  occur  in  May.  After  the  second  half  of  May,  the  mixed  layer  depth 
shallows  to  the  surface  and  remains  at  the  surface  over  the  summer.  After  the  second 
half  of  October,  mixed  layer  depth  again  deepens  until  end  of  the  year  to  120  m,  due  to 
vertical  overturning  caused  by  net  surface  cooling  plus  the  buoyancy  loss  attributable 
to  the  evaporation  rate. 
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Fig.   17.      The  Net  Heat  Flux  At  35  °  N  30  °  E 


turning  caused  by  net  surface  cooling  plus  the  buoyancy  loss  attributable  to  the  evapo- 
ration rate. 

Because  of  evaporation,  the  salinity  increases  at  fust  according  to  equation  (2.10). 
In  January  and  fall  the  predicted  mixed  layer  deepening  is  due  to  the  reduction  in  the 
downward  surface  buoyancy  flux,  causing  a  deeper  h.  As  evaporation  continues,  the 
surface  water  gets  saltier  and  denser  than  the  water  lying  immediately  below.  As  a  result, 
the  denser  surface  water  sinks  until  reaching  the  same  density.  Strong  positive  (E-P)  in- 
creases mixed  layer  salinity  by  increasing  the  numerator  of  the  equation  (2.10).  This  in- 
crease in  E-P  increases  the  mixed  layer  depth  by  decreasing  the  denominator  of  equation 
(2.39): 
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Fig.   18.      The  Mixed  Layer  Depth  With  Free  Convection  Only 
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See  Fig.  19  on  page  37  and  Fig.  20  on  page  38. 

In  the  central  Mediterranean,  the  mixed  layer  structure  is  very  similar  to  the  east- 
ern, with  only  slight  differences  at  corresponding  depths.  The  mixed  layer  deepens  to  197 
m  depth  at  this  location  in  January  and  it  stays  at  the  same  depth  between  February  and 
May  as  in  the  eastern  Mediterranean.  In  April,  the  mixed  layer  shallows  rapidly  to  the 
surface  and  then  stays  constant  until  the  middle  of  October.  After  October  the  mixed 
layer  depth  deepens  to  90  m  until  the  end  of  December.   See  Fig.  21  on  page  39. 
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Fig.   19.      The  Predicted  Salinity  With  Free  Convection  Only 


As  for  the  eastern  and  central  Mediterranean  mixed  layers,  the  western 
Mediterranean  mixed  layer  depth  is  well-correlated  with  the  net  heat  flux  variations. 
Going  from  the  eastern  to  the  western  Mediterranean,  the  rate  of  late  winter  and  early 
spring  mixed  layer  deepening  does  not  change  and  stays  constant  at  197  m  in  April. 
Conversely,  the  rate  of  mixed  layer  deepening  increases  at  the  end  of  the  year  from  60 
m  to  125  m,  as  we  go  from  eastern  to  western  Mediterranean.  See  Fig.  18  on  page  36, 
Fig.  21  on  page  39,  and  Fig.  22  on  page  40. 

B.      CASE  2  (THE  EFFECT  OF  THE  WIND  STRESS  ON  MLD  WITHOUT  E-P) 

Fig.  24  on  page  42  depicts  the  mixed  layer  depth  versus  time  in  the  eastern 
Mediterranean.    The  mixed  layer  depth  deepens  from  initial  depth  (95  m)  to  197  m  in 
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Fig.  20.      The  Predicted  Temperature  With  Free  Convection  Only 


January  in  the  forced  convection  plus  heat  flux  case  and  then  stays  constant  between 
February  and  May  at  197  m  as  in  the  free  convection  plus  heat  and  salinity  fluxes  case. 
In  April  there  is  a  sharp  shallowing  with  a  small  fluctuation.  Then  MLD  becomes 
constant  at  about  10  m  depth  until  September.  Between  October  and  January,  the  mixed 
layer  deepens  to  43  m.  In  April  there  is  a  large  sharp  decrease  and  then  a  small  sharp 
increase  in  mixed  layer  depth.  This  is  probably  due  to  a  rapid  changes  in  wind  stress  and 
surface  buoyancy  fluxes.  Also,  one  can  see  the  correlation  between  mixed  layer  temper- 
ature and  wind  velocity.  See  Fig.  25  on  page  43  and  Fig.  26  on  page  44.  From 
Fig.  27  on  page  45,  Fig.  24  on  page  42,  and  Fig.  28  on  page  46,  we  can  see  the  simi- 
larity between  the  mixed  layer  depth  plots. 
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Fig.  21.      The  Mixed  Layer  Depth  With  Free  Convection  Only 


In  the  central  Mediterranean  Sea,  the  mixed  layer  deepens  until  the  middle  of  Feb- 
ruary to  197  m.  Then  it  stays  constant  until  the  end  of  March.  After  March,  MLD 
shallows  sharply  with  a  small  fluctuation  in  April  until  July.  During  the  summer,  the 
mixed  layer  depth  is  less  than  6  m.  Due  to  the  stronger  wind  conditions  here  the  mixed 
layer  depth  deepens  between  October  and  January  to  56  m  depth. 

Mixed  layer  deepening  ends  in  February  in  the  western  Mediterranean  Sea  as  at  the 
other  locations,  due  to  increased  insolation.  Shallowing  starts  in  late  March,  and  it 
continues  until  July  at  around  10  m  depth.  Then  the  mixed  layer  remains  constant  until 
September  at  10  m  depth,  and  then  it  deepens  to  90  m  at  the  end  of  the  year  due  to  the 
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Fig.  22.      The  Mixed  Layer  Depth  With  Free  Convection  Only 


relatively  short  time  period  of  strong  winter  winds.   See  Fig.  34  on  page  52  and  Fig.  28 
on  page  46. 

C.      CASE  3  (COMBINED  FREE  AND  FORCED  CONVECTION) 

In  this  case  the  model  is  forced  by  the  combination  of  wind,  heat  fluxes  and  evap- 
oration. Between  October  and  March  the  mixed  layer  can  deepen  more  than  in  the  first 
two  cases  alone  due  to  the  combination  of  strong  wind  and  evaporative  conditions. 
During  the  fall  and  winter  cooling  period,  the  mixed  layer  cools  and  deepens.  Part  of 
the  cooling  is  due  to  surface  heat  loss,  and  part  is  attributable  to  the  entrainment  of 
underlying  cold  water.  The  seasonal  thermocline  weakens,  and  it  usually  vanishes  by 
March.  During  the  heating  period,  net  heat  gained  in  daytime  is  retained  through  the 
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Fig.  23.      The  Predicted  Temperature  With  Wind  And  Heat  Fluxes  Only 


night,  so  there  is  a  cumulative  gain.  Fig.  29  on  page  47,  Fig.  32  on  page  50,  Fig.  30  on 
page  48,  and  Fig.  31  on  page  49  show  mixed  layer  depth,  mixed  layer  temperature  and 
salinity,  and  incoming  solar  radiation  versus  time,  respectively,  for  the  eastern 
Mediterranean  Sea. 

D.     RESULTS 

The  most  significant  results  may  be  summarized  for  the  above  cases.  Whenever  the 
evaporation  becomes  largest,  the  mixed  layer  depth  becomes  deepest  in  the  case  of  free 
convection  forcing  only.  Depending  on  the  latent  heat  flux,  evaporation  is  largest  in 
winter  in  the  eastern  Mediterranean.  See  Fig.  33  on  page  51.  One  can  see  that  the  mixed 
layer  depth  reaches  its  maximum  value  after  the  late  winter  peak  in  the  latent  heat  flux 
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Fig.  24.      The  Mixed  Layer  Depth  With  Wind  And  Heat  Fluxes  Only 

and  therefore  in  evaporation.  There  is  a  good  correlation  between  latent  heat  flux  and 
mixed  layer  depth  at  every  station  in  the  Mediterranean  Sea. 

From  Fig.  26  on  page  44  and  Fig.  24,  the  wind  speed  becomes  largest  between 
November  and  March.  Mixed  layer  depth  also  reaches  its  deepest  value  of  197  m  at  this 
same  time  period  in  the  eastern  Mediterranean. 

Between  May  and  October,  mixed  layer  depths  are  at  a  minimum,  with  buoyant 
damping  associated  with  strongly  positive  Q0  exceeding  the  buoyant  production  due  to 
evaporation.   In  the  case  of  no  entrainment,  Wt  =  dhfdt  =  0,  and  B0  =  G0 , 


[<xg(^rCp)-pgS(E-PM  =  G0. 
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Fig.  25.      The  Predicted  Temperature  With  Wind  And  Heat  Fluxes  Only 


Therefore  the  second  term  on  the  left-hand  side  is  small  enough  to  make  the 
equation  (2.38)  positive.  The  deepening  in  the  mixed  layer  which  is  caused  by  wind 
forcing  is  almost  the  same  as  the  deepening  caused  by  the  evaporation  effect  between 
February  and  April.  However,  between  October  and  January,  the  deepening  in  the 
mixed  layer  caused  by  evaporation  exceeds  the  deepening  attributable  to  wind  stress 
alone  at  all  three  stations  in  the  Mediterranean  Sea.    See  Fig.   18  on  page  36,  Fig.  21 
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Fig.  26.      The  Wind  Speed  At  35  °  N  30  °  E 

on  page  39,  and  Fig.  22  on  page  40,  Fig.  24  on  page  42,  Fig.  27  on  page  45,  Fig.  28 
on  page  46. 

Over  the  year,  the  correlation  between  the  mixed  layer  depth  and  the  surface  heat 
flux  is  greater  than  between  mixed  layer  depth  and  the  evaporation  rate  or  between  the 
mixed  layer  depth  and  wind  stress.  In  particular,  the  net  heat  flux  minimum 
(November-April)  and  maximum  (May-October),  coincide  well  with  the  corresponding 
extremes  in  the  mixed  layer  depth. 

In  January-April  the  mixed  layer  depth  figures  show  clearly  the  domination  of 
evaporation  effects  over  the  wind  effect.  This  is  due  to  strong  evaporation  rate  over  the 
Mediterranean  Sea.    This  mixing  caused  bv  winds  is  both  free  and  forced.  The  resultant 
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Fig.  27.      The  Mixed  Layer  Depth  With  Wind  And  Heat  Fluxes  Only 

mixed  layer  depth  is  potentially  very  large,  and  can  be  a  cause  of  intermediate  and  deep 
water  formation. 

In  the  case  without  wind  mixing,  the  mixed  layer  depth  shallows  to  a  2  m  depth  in 
summer  time.  Thus  the  summer  wind  is  important  in  mixing  the  surface  layer,  and  can 
not  be  neglected.  This  is  due  to  the  domination  of  Q0  over  E-P  in  equation  (2.39).  In 
that  case  the  evaporation  process  becomes  less  insignificant  for  determining  h. 
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Fig.  28.      The  Mixed  Layer  Depth  With  Wind  And  Heat  Fluxes  Only 
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Fig.  29.      The  Mixed  Layer  Depth  With  Free  Convection  Plus  Forced  Convection 
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Fig.  30.      The  Predicted  Temperature  With  Free  Convection  Plus  Forced  Convection 
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Fig.  31.      The  Predicted  Salinity  With  Free  Convection  Plus  Forced  Convection 
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Fig.  32.      The  Incoming  Solar  Radiation  At  35  °  N  30  °  E 
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Fig.  33.      The  Latent  Heat  Flux  At  35  °  N  30  °  E 


Fig.  34.      The  Wind  Speed  At  40  °  N  05  °  E 


Fig.  35.      The  Net  Heat  Flux  At  40  °  N  05  °  E 
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Fig.  36.      The  Mixed  Layer  Depth  With  Free  Convection  Plus  Forced  Convection 


V.      CONCLUSIONS 

In  this  study  the  OPBL  model  was  employed  in  order  to  determine  the  effects  of 
seasonally  varying  variables  on  the  mixed  layer  dynamics.  The  surface  heat  flux,  wind 
stress,  and  evaporation  rate  are  provided  from  the  climatological  data.  For  a  variety  of 
cases,  the  model  is  solved  numerically  to  obtain  mixed  layer  depth,  mixed  layer  temper- 
ature and  salinity,  below- layer  temperature  and  salinity  profiles  for  a  whole  year.  Three 
numerical  experiments  are  conducted  to  demonstrate  the  relative  roles  of  wind  mixing, 
heat  flux  and  evaporation  in  determining  the  seasonal  evolution  of  the  mixed  layer  at 
three  sites  in  the  Mediterranean  Sea. 

In  the  case  of  free  convective  mixing  only  (first  case),  evaporation  and  surface 
cooling  processes  can  deepen  the  mixed  layer  to  197  m  in  January  in  the  eastern 
Mediterranean  Sea.  Free  convection  and  heat  flux  alone  can  deepen  the  mixed  layer  to 
same  depth  in  the  central  and  western  Mediterranean  as  in  the  eastern  Mediterranean. 
Over  the  summer  period,  free  convection  and  heat  flux  alone  can  not  explain  the  ob- 
served mixed  layer  depth  for  all  three  locations  in  the  Mediterranean  Sea.  In  the  fall  and 
early  winter,  free  convection  in  the  eastern  Mediterranean  Sea  can  deepen  the  mixed 
layer  to  60  m  depth.  In  the  central  and  western  Mediterranean  such  deepening  reaches 
90  m  and  125  m  depths  respectively.  Therefore,  in  general,  the  fall  mixed  layer  salinity- 
increases  and  the  temperature  decreases  due  to  excessive  evaporation  rate  over  the 
Mediterranean  Sea.  Then  in  winter  the  layer  entrains  the  less  saline  lower  layer,  bal- 
ancing the  salinity  increase  caused  by  evaporation.  The  positive  buoyancy  flux  caused 
by  increased  turbulent  kinetic  energy  due  to  strong  evaporation  starts  the  mixing  process 
in  fall  and  winter. 

When  the  wind  is  a  dominant  mechanism,  in  the  case  of  setting  net  evaporation  to 
zero  (second  case),  the  results  show  that  the  mixed  layer  deepening  reaches  197  m  in  late 
winter  at  all  three  stations. 

Therefore,  from  the  first  and  second  cases  in  late  winter  until  April,  the  wind  and 
evaporation  have  almost  the  same  effects  on  the  mixed  layer  deepening.  Ih  fall  and  early 
winter  until  January,  the  evaporation  process  is  dominant,  and  this  is  especially  true  in 
the  western  Mediterranean  Sea. 

This  experiment  showed  that  the  evaporation  is  weakest  in  summer,  and  is  most 
significant  in  late  fall  and  in  early  winter.  The  effect  of  wind  mixing  is  much  more  con- 


tinuous  in  time  and  space  than  is  the  effect  of  evaporation.  There  are  no  significant  wind 
effect  differences  between  the  three  stations. 

This  study  shows  that  for  the  seasonal  time  scale,  evaporation  is  an  important 
driving  mechanism  for  mixed  layer  deepening  in  the  Mediterranean  Sea. 

It  is  recommended  that  this  research  be  continued,  with  an  emphasis  on  the  forcing 
attributable  to  storms  and  shorter-period  fluctuations  in  the  surface  forcing. 
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